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ABSTRACT 

Aims. The aim of this work is to study the contribution of the Lyo- emitters to the star formation rate density (SFRD) of the Universe 
in the interval 2 < z < 6.6. 

Methods. We assembled a sample of 217 Lyo- emitters (LAE) from the Vimos-VUT Deep Survey (VVDS) with secure spectroscopic 
redshifts in the redshift range 2 < z < 6.62 and fluxes down to F~ 1.5 x lO^^^erg/s/cm^. 133 Lyo' emitters are serendipitous 
identifications in the 22 arcmin^ total slit area surveyed with the VVDS-Deep and the 3.3 arcmin^ from the VVDS Ultra-Deep survey, 
and 84 are targeted identifications in the 0.62 deg- surveyed with the VVDS-DEEP and 0. 16 deg^ from the Ultra-Deep survey. Among 
the serendipitous targets we estimate that 90% of the emission lines are most probably Lyo-, while the remaining 10% could be either 
[Oil] 3727 or LyQ-. We omputed the luminosity function and derived the star formation density from LAE at these redshifts. 
Results. The VVDS-LAE sample reaches faint line fluxes F(Lya) = 1.5 x lO^^^erg/s/cm^ (corresponding to L(Lya) ~ I0*^erg/s 
at z ~ 3) enabling the faint end slope of the luminosity function to be constrained to cc ~ -1.6 ± 0.12 at redshift z ~ 2.5 and to 

a 1.78*^012 at redshift ~ 4, placing on firm statistical grounds trends found in previous LAE studies, and indicating that sub-L, 

LAE (Liy-a :£ I0*~^erg/ s) contribute significantly to the SFRD. The projected number density and volume density of faint LAE in 
2 < z < 6.6 with F> 1.5 x 10""*erg/i/c/)i^ are 33 galaxies/arcmin^ and ~ 4 x lO^^Mpc"^, respectively. We find that the the observed 
luminosity function of LAE does not evolve from z=2 to z=6. This implies that, after correction for the redshift-dependent IGM 
absorption, the intrinsic LF must have evolved significantly over 3 Gyr The SFRD from LAE is found to be contributing about 20% 
of the SFRD at z = 2 - 3, while the LAE appear to be the dominant source of star formation producing ionizing photons in the early 
universe z ~> 5-6, becoming equivalent to that of Lyman Break galaxies. 

Key words. Cosmology: observations - Galaxies: fundamental parameters - Galaxies: evolution - Galaxies: formation 



1. Introduction 

The Lye line is the strongest hydrogen emission line in the 
Universe, and it is observed in the optical range for galaxies 
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at z > 2. It has thus naturally been used to search for high- 
z galaxies (Partridge&Peebles 1967; Djorgovski et al. 1985; 
Cowie &Hu 1998). 

The Lya emission in galaxies is thought to be produced by 
star formation, as the AGN contribution to the Lyor population 
at z < 4 is found to be less than 5% (Gawiser et al. 2006; 
Ouchi et al. 2008; Nilsson et al. 2009). However, the physical 
interpretation of the observed Lya flux is not simple, because 
the Lya photons are resonantly scattered by neutral hydrogen. 
Lya photons can therefore be more attenuated than other UV 
photons, and they have an escape fraction that can depend on 
the spatial distribution of neutral and ionized gas, as well as 
on the velocity field of the neutral gas (Giavalisco et al. 1996; 
Kunth et al. 1998; Mas-Hesse et al. 2003; Deharveng et al. 2008 
and references therein). 
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hya emitters (LAE) observed up to now are forming stars at 
rates of ~ 1 lOM^yr-^ (Cowie&Hu 1998; Gawiser et al. 2006; 
Pirzkal et al. 2007), and they have stellar masses as low 
as 10^ -;- lO^Mo and ages < 50M>t (Pirzkal et al. 2007; 
; Gawiser et al. 2007; Nilsson et al. 2009). However, 
Nilsson et al. (2007) find ages between 0. 1 and 0.9 Gyr, and more 
recently Pentericci et al. (2007) and Finkelstein et al. (2009) 
point out that Lya galaxies are a more heterogeneous family than 
young star-forming galaxies: they also find Lya emitters with old 
stellar populations (ages of ~ 1 Gyr) and a wide range of stellar 
masses (up to IO'^Mq). 

Interestingly, a class of Lya galaxies with rest-frame equiva- 
lent width EW > 240A has been found (Malhotra&Roads 2002; 
Shimasaku et al. 2006). Galaxies with such large EW cannot 
be explained by star formation with a Salpeter IMF, but must 
have a top heavy IMF, a very young age < 10^ yr and/or a very 
low metallicity. Many of these large EW objects are spatially 
extended, and thus are good candidates to be cooling clouds or 
primeval galaxies (Yang et al. 2006; Schaerer 2002) and can give 
interesting clues about the first stages of star formation. 

Together with studying the properties of Lya galaxies at 
different redshifts, it is important to study the evolution of 
their luminosity function, comparing large and complete sam- 
ples of Lya galaxies at diff'erent redshifts. The most com- 
mon technique used so far has been to build large samples 
from imaging in narrow band filters tuned to detect Lyor emis- 
sion at z ~ 2 ^ 9 (Hu et al. 2004; Cuby et al. 2003; 
Tapken et al. 2006; Kashikawa et al. 2006; Gronwall et al. 2007; 
Murayama et al. 2007; Ouchi et al. 2008; Nilsson et al. 2009; 
Guaita et al. 2010). Blank field spectroscopy has been used 
blindly to search for Lyo' emitters in deep HST-ACS slit- 
less spectroscopic observations (MaUiotra et al., 2005) or 
slit spectroscopy (van Breukelen, Jarvis & Venemans 2005; 
Martin et al. 2008; Rauch et al. 2008; Sawicki et al. 2008), with 
Ranch et al. (2008) exploring the faintest emitters. 

The general consensus today is that the apparent 
luminosity function of Lya galaxies, that is the non- 
IGM corrected luminosity function, does not evolve at 
z ~ 3 6 (Rhoads&Malhotra 2001; Ouchi et al. 2003; 
van Breukelen, Jarvis & Venemans 2005; Shimasaku et al. 2006; 
Murayama et al. 2007; Gronwall et al. 2007; Ouchi et al. 2008; 
Grove et al. 2009). However, this conclusion is drawn 
from small samples, as the narrow band imaging tech- 
niques sample only thin slices in redshift (typically Az=0.1) 
and needs to be spectroscopically confirmed. Extensive 
spectroscopic follow-ups of narrow band Lya candidates 
have been carried out in recent years, gathering hun- 
dreds of spectroscopic confirmations between z ~ 2 and 
z ~ 7. However, the spectroscopic coverage rarely reaches 
30 -;- 50% of the photometric sample (Murayama et al. 2007; 
Gronwall et al. 2007; Ouchi et al. 2008), but usually is 
much lower (Kashikawa et al. 2006; Nilsson et al. 2007; 
Matsuda et al. 2005). Moreover, current blind spectroscopic 
surveys sample only small areas to relatively shallow fluxes. 
It is likely that the apparent lack of evolution is a coincidence 
of the evolving intrinsic Lya LF combined with an evolution 
of the intergalactic medium absorption with redshift (e.g. 
Ouchi et al. 2008). Firm conclusions about the possible evo- 
lution of the luminosity function have not yet been secured. 
Moreover, existing spectroscopic and narrow band samples are 
not sufficiently deep to constrain, even at intermediate redshift 
(z ~ 3), the slope of the luminosity function. 

Measuring the luminosity function in turn enables to com- 
pute the luminosity density and star formation rate density evo- 



lution under a set of well constrained hypotheses. The contri- 
bution of Lya to the total star formation rate is yet not robustly 
measured mainly because the faint end slope of the luminosity 
function remains poorly constrained. 

In this paper, we present the results from a very deep blind 
spectroscopic survey search of Lya emitters over an unprece- 
dented large sky area. We looked for the serendipitous detection 
of Lya emission in the slits of the VIMOS VLT Deep Survey, 
concentrating on the VVDS-Deep and VVDS-Ultra-Deep sur- 
veys reaching up to 18h of integration on the VLT- VIMOS. We 
describe the spectroscopic and photometric data and the associ- 
ated selection function in Section 2. The search for Lya emitters 
and the final sample are presented in Section 3, and we discuss 
its properties in Section 4. The luminosity function calculation 
is discussed in Section 5, and the star formation rate density is 
derived. We discuss these results and give a summary in Section 
6. 

Throughout the paper, we use and AB magnitudes and a stan- 
dard Cosmology with £Im = 0.3, Qa = 0.7 and h = 0.7. 

2. Search for serendipitous emission lines in the 
VVDS Deep and Ultra-Deep spectroscopic 
surveys 

The Vimos-VLT Deep Survey (VVDS) exploited the high mul- 
tiplex capabilities of the VIMOS instrument on the ESO-VLT 
(Le Fevre et al., 2003) to collect more than 45000 spectra of 
galaxies between z ~0 and z ~ 5 (Le Fevre et al., 2005; 
Garilli et al. 2008). In the VVDS-Deep 0216-04 field, more flian 
~ 10000 spectra have been collected for galaxies with Iab ^ 24, 
observed with the LR-Red grism across the wavelength range 
5500 < A < 9350A, with integration times of 16000 seconds. 
In addition, the VVDS Ultra-Deep (Le Fevre et al. 2010, in 
preparation) collected ~1200 spectra for galaxies with Iab ^ 
24.75, obtained with LR-blue and LR-red grisms, with integra- 
tion times of 65000 seconds for each grism. This produces spec- 
tra with a wavelength range 3600 < A < 9350A. For both the 
Deep and Ultra-Deep surveys, the slits have been designed to be 
1" in width, providing a good sampling of the 1" typical seeing 
of Paranal, and between ~4" and ~15" in length, allowing good 
sky determination on both sides of the main target. The resulting 
spectral resolution is R =^ 230 - 250 for both the LR-blue and 
LR-red grisms. The observations of the VVDS-Deep and Ultra- 
Deep have been taken in different observing runs split typically 
between 3 to 5 nights respectively. Given the Paranal seeing vari- 
ations, the seeing quahty of the final observations spans from 
0.5" to 1.2"FWHM. 

In this work, we take advantage of the fact that the VIMOS 
spectrograph produces a spectrum for the whole piece of sky 
covered by each sht. If a Lya emitting galaxy with a redshift 
compatible with our wavelength range serendipitously falls in 
the slit, the Lya line will appear in the 2-d spectrum. Obviously 
other line emitting galaxies at redshifts such that one or more 
of their emission line spectrum falls in the observed wavelength 
range wiU also be detected. 

During the processing of the 2-d spectra of the VVDS main 
targets we identified that a population of serendipitous emission 
line objects was present. We therefore subsequently performed 
a systematic search for serendipitous emission line galaxies in 
the ~ 8000 Deep and ~ 1200 Ultra-Deep spectra. This was then 
followed by unambiguous redshift identification of all the line 
emitting galaxies as described in Section 3.2. A significant num- 
ber of VVDS primary targets are also Lya emitting galaxies (Le 
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Fevre et al., 2010, in preparation), we added them to build a more 
complete LAE sample. 

2.1. The dataset 

The dataset consists of three sets of VIMOS pointings: the Deep, 
the Ultra-Deep blue, and the Ultra-Deep red. Each pointing con- 
sists of 4 separate quadrants, each containing approximately 100 
slits. We summarize each dataset below: 

- Deep dataset: it consists of 20 pointings, for a total of 80 
quadrants and about 8000 slits; the main VVDS targets have 
been selected from an area of 0.62 deg- to have 17.5< 
Iab <24. The spectra have been obtained with the LR-Red 
grism, providing a wavelength range 5500 < A < 9350A. 
The exposure times are about 16000 seconds, allowing to 
reach fluxes as low as ~ 5 x 10^^^ e rg / cm^ /s (see details 
below). The total sky area serendipitously covered is 22 
arcmin^. 

- Ultra-Deep dataset: it consists of 3 pointings, for a total of 
12 quadrants and about 1200 slits; the main targets have 
been selected to have 23 < Iab <24.75 from an area of 
0.16 deg^. Each slit placed on a primary VVDS target has 
been observed twice, once with the LR-Blue grism, and 
once with the LR-Red. For the blue spectra, the spectral 
range is 3600 < A < 6800A, and for the red ones it is 
5500 < A < 9350A. Thus, for each of the 1200 slits, the 
Ultra-Deep blue and red spectra overlap for about 1300A. 
We did not try to combine spectra in the overlapping regions, 
because blue and red spectra have been observed under dif- 
ferent sky background and atmospheric seeing conditions. 
As we will show in the next Sections, 13 Lya emissions 
have been found in the overlapping regions. As we handled 
the blue and red spectra separately this allows us to com- 
pare a posteriori the detection rate, as well as redshifts and 
fluxes. This is described in the next sections, respectively re- 
ferring to the two samples as 'Ultra-Deep blue' and 'Ultra- 
Deep red'. For each of the blue and red spectra the expo- 
sure time is 65000 seconds, allowing to reach a flux limit 
of ~ l.SxlO^^^erg/cm^/s (details below). The total sky area 
serendipitously covered is 3.3 arcmin^. 

The noise of VIMOS spectra varies as a function of wave- 
length, as a result of a combination of different effects. First, the 
combination of the efliciency of VIMOS optics and the CCD 
quantum efficiency depends on the wavelength of the incident 
light. Second, the OH airglow sky emission lines are in a suite 
of bands which are dominating the background for A > 7500A, 
and are not easy to subtract from the combined object-Hsky spec- 
tra in low to medium resolution spectra. As a result, spectra 
are noisier at the wavelengths of the OH bands. An added un- 
certainty on the flux level is also sometimes present at the po- 
sition of serendipitous Lyo' emission because the background 
subtraction performed by the VIMOS data processing pipeline 
VIPGI (Scodeggio et al., 2005) is optimized using the a-priori 
knowledge of the position of the VVDS primary target. As the 
serendipitous target is offset from the main target, the low or- 
der polynomial fit used to remove the background is made using 
some background points including the emission line, therefore 
artificially lowering the line flux. Third, the thin E2V detectors 
on VIMOS produce fringing which strongly affects the part of 
the spectra at A > 8000A. As the VVDS has been using a low 
resolution R ~ 230 to maximize the number of objects in the 
survey, this produces a blend of the OH sky emission features. 
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Fig. 1. Detection limits for the hya line, as a function of wavelength 
and redshift, for the three subsamples: Deep (green dashed line), Ultra- 
Deep red (red dotted line) and Ultra-Deep blue (blue line). The shaded 
regions indicate part of the spectra dominated by bright OH sky lines, 
where the rms of the background becomes too high to enable a secure 
line detection. 



globally increasing the background noise at the position of the 
OH bands. As a result of these eff'ects, the flux limit varies as 
a function of wavelength, hence of redshift. We determined the 
theoretical flux limit, separately for 'Ultra-Deep blue', 'Ultra- 
Deep red' and Deep, by measuring the typical rms at different 
wavelengths in the 2 - d spectra. Empirically, we decided to set 
the limit to 3cr in at least 5 contiguous pixels. The limits for 
the three subsamples are shown as a function of wavelength in 
Figure 1 . It can be noted that the background, even for each of 
the three subset, is changing as a function of the wavelength and 
redshift. As a result of the bright OH skylines at A >7500A, a 
large part of the redshift range at z > 5.5 is in practice inaccessi- 
ble. However, there are at least two very clean windows around 
z-5.1 and z=6.5, with a very low background coiTesponding to 
the absence of OH emission. The faintest flux limits can be ob- 
served in the range 4200 < A < 7600A, reaching down to very 
faint flux levels of ~ 1.5xl0"'^erg/cm^/5 at 3cr. The combi- 
nation of sky area covered, wavelength range, and depth is un- 
precedented. 

2.2. Photometry 

The VVDS 0216-04 field benefits from extensive deep pho- 
tometry. The field was first observed with the CFH12K cam- 
era in BVRI bands (Le Fevre et al., 2004; McCracken et al., 
2003), in the U band (Radovich et al. 2004) and K band 
(lovino et al. 2005). More recent and significantly deeper obser- 
vations have been obtained as part of the CFHT Legacy Survey 
in ugri and z' (Goranova et al. 2009; Coupon et al. 2009), and as 
part of the WIRDS survey in J, H and K bands. The CFHTLS 
observations reach 5cr point source limiting magnitudes in i band 
of lAB - 28 while the WIRDS observations reach a 3cr point 
source limit of Ksab - 23.5 (McCracken et al., in preparation). 
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3. The Lya sample 

3.1. Serendipitous emission lines identification 

The final database consists of ~ 1200 2-d spectra observed with 
both red and blue grism (Ultra-Deep dataset), and ~ 8000 2-d 
spectra observed with the red grism. For each of the 92 quadrants 
observed, the reduction code (VIPGI, Scodeggio et al. 2005) 
produces a frame in which the ~100 2-d spectra are arranged 
one above the other, with the spectral direction along the x-axis. 
With the purpose of identifying serendipitous lines, we built a 
data processing procedure that automatically masked out all the 
continua and the bad features (bright residuals and OH skylines). 
The unmasked area is then used to search for emission Unes, and 
to compute the effective area covered by each slit on the sky. An 
object falling serendipitously in a slit will produce a spectrum 
with a combination of continuum and emission lines. For faint 
objects with strong lines, the continuum is not detected, and we 
are left with emission lines with a Une profile as produced by the 
spectrograph. For point source objects, VIMOS produces a Une 
profile along the slit (in the spatial direction) which is almost 
identical to the seeing profile (FWHM~ 1" or about 5 pixels). 
This is the projected equivalent slit size resulting from the ob- 
ject size convolved with the seeing profile and the slit profile 
on the spectral dispersion direction for a maximum of about 5 
pixels corresponding to the 1 arcsecond slit width projection (or 
about 30A spectrally). The 2D Une profiles produced are thus 
similar to point source images. For compact objects smaller than 
the atmospheric seeing disc, and if the seeing is smaller than the 
1 arcsec slit width, the observed line profile in the spectral di- 
mension will be dominated by the projected seeing profile. For 
extended objects the Une profile becomes wider in the spatial 
direction leading to oval shaped emission Unes. As these spec- 
troscopic line images are similar to those of stars or galaxies 
on deep images, we blindly run SExtractor on the masked 2D 
spectra to automatically find all emission lines. This was fol- 
lowed by a visual inspection by two authors (PC and OLF) of 
all the objects detected by SExtractor, each producing a catalog 
of candidates. A further inspection has been conducted in order 
to find possible lines missed by SExtractor: we noted that lines 
too close to a bright continuum were not found by SExtractor, 
probably because of deblending problems. These two indepen- 
dent lists were then jointly examined by these two authors to 
decide if a candidate was to be kept, with more than 85% of the 
Usts being unambiguously and independently selected by both 
authors. A fiducial catalog of bona-fide emission lines has then 
been built from this comparison. A more careful determination 
of the incompleteness and contamination will be presented in the 
next sections. In particular, the completeness is assessed using 
simulations in Section 3.4, and an estimate of the contamination 
is given in Section 3.2. 

We obtain a total sample of 133 serendipitous emission line 
emitters, 105 found in the Ultra-Deep survey and 28 in the Deep. 
Most of these galaxies have no evident continuum appearing in 
the spectra. We show a random pick of 8 serendipitous Unes in 
Figure 2 and 3.1. For each line, we show both the 2-d and 1-d 
spectra. 

3.2. Lya identification 

Once the fiducial emission Une catalog is built, we have to unam- 
biguously assign a rest-frame wavelength to any given emission 
line. If one emission line is identified in a contiguous spectrum 
covering 3600 to 9350 A, there are several possibiUties for the 



line identification assuming a star-forming galaxy spectrum: be- 
low 3727A the line is most probably Lya , above 3727A the 
emission Une can be Lya , [OII]3727A, H/3, [OIIIJ5007A (as- 
suming [OIII]4959A is not detected, [OIII]5007/[OIII]4959^3), 
or Ha. In the spectrum of normal star forming galaxies, [Oil], 
Ufi and [OIII] are most often present. Therefore, as a first ap- 
proach we assumed that the emission lines could be [Oil] at 
3727 A, HyS at 4861 A, [OIII] at 5007 Aor Ha at 6563 A, and 
we checked for other expected emission Unes at different wave- 
lengths in our spectral window. Using the typical Une ratios for 
star-forming galaxies of average metallicity we estimate the flux 
of the other expected lines, and check whether they could be de- 
tected above our flux limit. However, the line flux ratios strongly 
depend on the galaxy metallicity. For example, typical ranges 
of line ratios are: [OIII]/[OII]= 0.1 ^ 5, Ha/[OII]= 0.3 - 5 and 
H/3/Ha = 0.1 H- 0.8 (Lamareille et al. 2006; Maier et al. 2006; 
Tresse et al. 2007; Cowie&Barger 2008, Kewley&ElUson 2008). 
So, we decided to use the foUowing Une ratios, typical of galax- 
ies wifli average metalUcity: [Oni]/[On]~0.35,Ha/[On]~l and 
HyS/Ha ~0.35. Cases of extreme metallicity are discussed below. 

Narrow-line type-II AGN with classical line ratios would 
be straightforward to identify given the broad wavelength range 
of our spectra. If we were to assign the single emissions Unes 
to lines like CIV-1549, CIII-1909, or MgII-2800 we would ex- 
pect to see other emission lines in our observed range either 
towards the UV (Lya, CIV, CIII) or towards the red (Mgll, 
[Oil], [Nelll]). Broad line type-I AGN are also relatively easy 
to identify at our spectral resolution. We have not found any 
type-I or type-II AGN in our data, consistent with expected 
AGN counts in the area sampled (e.g. Bongiomo et al., 2007; 
Lamareille et al., 2009). 

As the spectral window spans from 3600A to 9350A for the 
Ultra-Deep subsample and from 5500A to 9350A for the deep 
subsample, we analyze these samples separately. 

- Ultra-Deep: 2 single lines in the Ultra-Deep are at 
A <3727A, and then can likely be identified as Lya , while 
103 are at A >3727A. Among the latter, 93 Unes are at 
A <6920A. If they were [Oil], we would to have [OIII] and 
Hy6 in our spectral window. We estimate, using the line ratios 
reported above, that 75 of these 93 lines should be accom- 
panied by either detectable [OIII] or detectable Hfi or both, 
but as they are not detected they are therefore most proba- 
bly to be identified with Lya . For the 10 emission lines at 
A >6920A, this argument cannot be used, because the possi- 
ble [OIII] and Ha emission would Ue out from our spectral 
window. Using similar arguments, we can infer that none 
of the Ultra-Deep lines are H/?, [OIII] or Ha, because they 
would be always accompanied by at least one other line in 
our spectral window with a flux brighter than the flux Umit. 

- Deep: The deep subsample contains 28 serendipitous Unes. 
Again, all of them can be in principle interpreted as both Lya 
or [Oil] . Since for the Deep the spectral window is smaller 
than for the Ultra-Deep, the same argument as for the ultra- 
deep tells us that only 7 of these lines are not [Oil]. For the 
other 2 1 , either a possible [OIII] emission is outside the spec- 
tral range, or it would be too faint to be detected. However, 
we can conclude also that the bulk of them are not Ha, [Olll] 
or H/3, because at least one other detectable Une would be de- 
tected in our spectra. 

From this first analysis, we can conclude that basically none 
of the serendipitous lines in the sample are Ha, H^ or [OIII], 
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Fig. 2. Random selection of LAE candidates from the Ultra-deep blue dataset, ordered by increasing redshift. For each of them we show the full 
2-d (top) and 1-d (bottom) spectra. The position of Lyo- is indicated by an arrow on both the 2-d and 1-d spectra, and by a label reporting the 
redshift on the 1-d spectrum. The observed wavelength scale at the bottom of the 1-d spectrum is transformed into the redshift of Lyo- at each 
wavelength on the top. 
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Fig. 3. Same as figure 2, for objects detected in thie Ultra-deep red and in the Deep datasets. 
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Fig. 4. top panel: For tiie full sample, the distribution of rest-frame EW 
if the lines are identified as Lya (black histogram) or as [OII]3727 (red 
histogram). The filled histograms show objects without a detected con- 
tinuum, for which the EW is a lower limit, lower panel: the same as 
the top panel, but for the 49 lines identified as ambiguous. The blue and 
green histograms represent respectively the rest-frame EW if the lines 
are identified as Lya or as [Oil] 3727. 



that 84 among the 133 bona fide serendipitous lines are not 
[Oil], and that most of the line identification uncertainty for 
the other 49 emission lines is between [Oil] and Lya . We 
note that these findings will not change too much if we had 
used much lower values for the metallicity. Assuming, for ex- 
ample, [OIII]/[OII]~10, if the serendipitous line was assigned to 
[Oil], we should always have detected [OIII] for all lines with 
A < 6920. 

We can perform a further check in order to evaluate the con- 
tribution of Lye and [Oil] emission to the observed population. 
In particular, we measured the observed equivalent width of the 
lines to produce the distribution of the rest-frame EW for each 
of the two cases, as presented in Figure 4, for all the 217 lines 
as well as for the 49 ambiguous lines. While the distribution of 
the rest-frame EW for the lines if they are Lya is mostly below 
150A, we find that if the lines were assigned to [Oil], ~70% 
of the distribution would be with EW[OII]> lOOA. Among the 
49 galaxies that cannot be unambiguously classified as either 
Lyo- or [Oil], 37 would have EW[OII]> lOOA. We compare this 
distribution to the distribution of EW[OII] observed for galax- 
ies with Mb - -18 in the VVDS-Deep (Vergani et al. 2008) 
in the same redshift range 0.86 < z < 1.5, and we note that 
90% of the normal galaxy population has EW[OII]< lOOA. 
At z - 1, galaxies with EW[OII]> lOOA have only rarely 
been observed in deep spectroscopy surveys, and none with 
EW> 150A (Hammer et al., 1997; Vergani et al., 2008). Since 
we expect to have just 5 lines with EW[OII]> lOOA (10% of 
49), we can conclude that for about 32 (37-5) the identification 
is very likely to be Lya . So, at the most we can identify the 
line as [Oil] for 17 galaxies (49-32). We therefore conclude that 



Fig. 5. Combined spectra for all serendipitous LAE in 4 redshift do- 
mains, centered on Lya. The asymmetry of the line profile is clearly 
seen in each redshift range when comparing to a gaussian fit of the red 
wing of the line (long-dashed line). For each combined spectrum we 
report also the skewness calculated in the range 1201 < X < 1231 A. 
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Fig. 6. Combined spectra for 5.9 < z < 6.62. The top panel shows the 
continuum bluewards of the Lyor line. The bottom one focuses on the 
line, showing the asymmetric profile. The value for the skewness of the 
line is also reported. 



about 10% of the single emission lines could be assigned to ei- 
ther [Oil] or Lya, while more 90% of them are most likely Lya. 

Another check on the line identification with Lya is to look 
at the line profile. Lya line profiles are most often asymmetric, 
with a truncation towards the blue because of absortion by the 
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intergalactic medium, and an extended red wing because of com- 
plex velocity fields (Shapley et al., 2003; Ouchi et al., 2008). On 
individual objects, the line detection is not always of sufficient 
S/N to detect a line asymmetry, instead we produced the average 
normalized Lya line profiles for LAE in different redshift ranges 
up to z ~ 6, as shown in Figure 5 (the description of the combi- 
nation technique and a more detailed analysis of the combined 
spectra are presented in Section 4.3). The asymmetry of the pro- 
file is clearly seen at all redshifts, even in the lower redshift bins. 
For each combined spectrum we measured also the skewness of 
the line in the range 1201 < A < 1231 A , that measures the 
degree of asymmetry of a given distribution. If the line profile 
would be perfectly described by a gaussian, that is by definition 
symmetric, the skewness would be zero. On the other hand, a 
positive value for the skewness indicates a distribution skewed 
towards red wavelengths. As Figure 5 shows, in each redshift 
bin between 2 < z < 5.9 the skewness is positive and is increas- 
ing with redshift as expected from IGM absorption, reinforcing 
the primary identification with Lya . 

There are 6 galaxies with a single enoission line in the atmo- 
spheric windows 8450 <A< 8620 and 9000 <A< 9300 clean 
of OH sky emission lines, and they observed equivalent widths 
540 < EW < 1250A. If ffiese fines were [OIIJ3727 or Ha, ffie 
rest EW would be 220 < EW < 510 or 400 < EW < 880A with 
a mean EW=3 15 or 550A respectively, clearly outside the range 
of EW for I Oil] or Ha emitting galaxies at z ~ 1 . No sign of the 
broad component of an AGN which could produce high Ha EW 
has been observed. We also produced the fine profile of the 6 
galaxies in this redshift range, reported in Figure 6. The resulting 
profile is clearly asymetric, as confirmed by the positive value 
for the skewness. As other single emission line candidates with 
strong observed EW are even more unlikely than [Oil] 3727 or 
Ha and should show other emission Unes in our observed wave- 
length domain, we argue that the most likely possibiUty is that 
these emission lines are Lya with 5.96 < z < 6.62, making them 
some of the most distant galaxies found to date. 

In summary, among the 133 emission lines identified 
serendipitously in the 2D spectra, we conclude that 124 are most 
likely Lya , while a low fraction of 14 (10%) can be either Lya 
or [Oil] . In the following, we assume that all galaxies with am- 
biguous line assignment are Lya . 

3.3. Final sample 

We completed the LAE sample by adding the primary VVDS 
spectroscopic targets that unambiguously have Lya in emission. 
We found 70 Lya emitters among the ~ 1200 Ultra-Deep pri- 
mary targets, all at redshift z < 3.5, and 14 among the ~8000 
Deep targets. The Deep survey is actually limited to A >5500A, 
hence it could not see Lya at z < 3.5. All these galaxies have a 
high quality flag (larger than or equal to 2; Le Fevre et al., 2010, 
in prep.) for the redshift, since they have other spectral features 
that allow to unambiguously identify the redshift and therefore 
classify them as Lya emitters. They show in fact OI at 1303 A , 
CII at 1333A , SilV at 1397A , CIV at 1549A, and sometimes 
cm at 1909 A . The final sample is flierefore made of 217 LAE, 
including 133 serendipitous LAE, and 84 LAE from the primary 
VVDS targets. The sample is summarized in Table 1 . 

For each emission line, we carefully measured the position 
of the line (and thus the redshift), the total flux in the line and 
the continuum. The continuum, in units of F^, is measured as 
close as possible to the red wing of the Lya line, by averaging 
the counts between 1230 and 1250A. From these fundamental 





U-DEEP 


U-DEEP 


DEEP 


total 




blue 


red 






Serend. 


94' (46) 


24' (7) 


28(13) 


133(66) 


Targets 


70 





14 


84 



Table 1. The final sample of Lya emitters, divided into serendipitous 
objects and spectroscopic targets. The number in parenthesis for the 
serendipitous galaxies indicates the objects with optical counterpart. 
Note:' 13 lines have been observed in both U-DEEP blue and U-SEEP 
red 



quantities, the observed and rest-frame equivalent widths and the 
luminosities are derived: 

L = AnFdiizf EWrest = F/C-^ (1) 

1 +z 

where L and F are respectively the luminosity and the flux in the 
line, z is the redshift of the line, di{z) is the luminosity distance 
at the redshift of the line, EWrest is the rest-frame Equivalent 
Width and C is the continuum around the fine. Obviously, when 
the continuum is not detected, the derived equivalent widths are 
just a lower limit, and we used 1 sigma of the background as the 
continuum value. 

Moreover, each luminosity can be converted to a star forma- 
tion rate (SFR) combining the classical prescription to derive the 
SFR from Ha Luminosity by Kennicutt (1998) and the expected 
ratio between Lya and Ha ffux. In particular: 

SFRiMeyr-') = Lnaiergs s-')/il.26 x 10^^), (2) 

with Ha flux corrected for internal extinction and valid for a 
Salpeter IMF. The conversion between Ha and Lya luminosity 
has been theoretically derived by Brocklehurst 1971 (assuming 
case B recombination): 

Llya = ^-TLna, (3) 

and it does not account for dust and escape fraction corrections. 
Combining the two gives: 

SFRiMeyr-^) = Uyjergs s-')/(L\ x lO^^). (4) 

In Figure 7 we report the ffux and the position (both wave- 
length and redshift if the Unes are identiffed as Lya ) for the 
three surveys separately. It can be noted that the Ultra-Deep blue 
is providing the biggest sample with 164 emission lines galaxies 
(94 serendipitous lines and 70 targets of the spectroscopy), while 
the Ultra-Deep red and the Deep contribute with respectively 24 
and 42 emission Unes galaxies. Note also that in a few instances 
we measure line fluxes that are weaker than the nominal flux 
limit at that wavelength: this is because the background, at each 
wavelength, is measured as the average of the backgrounds in 
the different quadrants and pointings. For some pointings, the 
background is lower than the average, and thus fainter lines can 
be detected. 

3.4. Completeness simulations 

As we saw in Fig. 1, the background of our spectra has a com- 
plicated structure, as there are strong variations with wavelength. 
To perform a statistical analysis of the number density evolution 
of Lya galaxies, it is important to carefully determine, for each 
wavelength (and thus for each redshift), what is our capabiUty to 
recover a line with a given luminosity. 
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Fig. 7. For each sub-dataset, we show the flux versus redshift together 
with the flux limit as a function of the redshift and wavelength (solid 
lines). The flux hmit is empirically measured on the 2-d spectra. Filled 
and open circles represent respectively serendipitous Lya galaxies and 
targets with Lya emission. 



We performed a Monte Carlo simulation, building a catalog 
of 1000 fake lines for each dataset, for a total of 3000. Each line 
has a known flux, randomly extracted from a flat distribution be- 
tween Fiya) = 1 X 10"'** and Fiya) = 1 x 10^^''. The redshift is 
randomly extracted from a flat distribution between z=2 and and 
z=4.7 (for blue spectra) and between z=3.7 and z=6.7 (for red 
spectra). The fake lines are then added at the wavelength cor- 
responding to the known redshift, in a randomly extracted slit. 
Since emission lines in our sample are not or barely resolved 
either in the spatial or in the spectral directions (see next sec- 
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Fig. 8. Completeness in line identification from Monte Carlo simulation 
as a function of the flux, for the three subsamples: Deep (solid line), 
Ultra-Deep red (dotted line) and Ultra-Deep blue (dashed line). 



tions), we chose to use simple gaussian profiles with FWHM of 
l"x30A. 

We attempted to recover these simulated lines, without a- 
priori knowledge of redshifts and fluxes, by applying the same 
procedure used on the real 2-d spectra, described in Section 3.1. 
The results are reported in Figure 8 presenting the completeness 
as a function of the line luminosity and of the redshift, for the 
three subsamples. We can see that Ultra-Deep blue & red and 
Deep are 50% complete respectively at F^ya ~ 2 x 10"^^ and 
-8x10" 



-18 



3.5. Slit flux loss 

Since the slits used have sizes comparable to the seeing of the 
images, part of the line flux cannot be recovered and will be al- 
ways lost. The flux loss depends on the offset between the galaxy 
and the slit center positions, on the size of the galaxy with re- 
spect to the slit width and on the atmospheric seeing. For an ob- 
ject perfectly centered in the slit, the flux loss will be minimal, 
but not zero. If the position of the object in the slit is known, as in 
the case of the objects with an optical counterpart, it is relatively 
easy to compute the flux loss. However, a significant number of 
galaxies in our sample do not have an optical counterpart within 
+0.5" from the expected position of the object along the slit and 
within ±2" off axis in the direction perpendicular to the slit (see 
Section 4.2). 

As we do not know for these galaxies their positions inside 
the slits, we cannot compute object by object the actual flux loss. 
This problem is common to all surveys of serendipitous Lya 
emitters (e.g. Rauch et al. 2008; Lemaux et al. 2009), and is 
usually solved with a statistical approach for the objects with no 
detected counterpart. In this Section, we derive the slit flux loss 
as a function of the position of the object in the slit, so that we 
can correct the line flux for those objects with an optical counter- 
part, and then we estimate a statistical correction for the objects 
without any optical counterpart. 

We start from the consideration that the bulk of our Lya 
galaxies are unresolved in the spatial direction (see Section 4.2), 
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Fig. 9. Histogram of the percentage of recovered flux, for a set of a 
3000 simulated lines randomly offset with respect to the center of the 
slit (solid line) and for the serendipitous lines (dotted line). The vertical 
lines show the median for the two distributions, around 55% of the re- 
covered flux. The small inset shows the recovered flux as a function of 
the offsets for the simulated lines, for the 0.7, 0.9 and 1.2 arcsec FWHM 
lines (respectively in red, black and blue). 



having sizes comparable with the seeing of the observations. The 
FWHM distribution in the spatial direction of the lines in our 
sample spans from 0.5 to 1 .2 arcseconds, with 90% of them hav- 
ing FWHM=0.7-h1.2 arcsec (see Figure 11). This implies that 
we can model them as gaussians, with fixed FWHM. We built a 
simulation in which an emission line object with a given FWHM 
is placed in a 1 arcsec slit at diff'erent positions with respect to the 
center of the slit until it is completely outside of the slit width. 
The small inset in Figure 9 shows the fraction of recovered flux 
as a function of the off'set between the objects and the slit, for 
three seeing conditions FWHM= 0.7, 0.9 and 1.2 arcsec (cor- 
responding to cr ~0.3, 0.38 and 0.5 arcsec). If the galaxy was 
perfectly centered in the slit, only ~ 70, 80 and 90% of the flux 
would be recovered, respectively for the 1.2, 0.9 and 0.7 arcsec 
FWHM lines; on the other hand, if the galaxy would be placed 
at 2cr=0.8 arcsec outside the slit, only 20% of the light would be 
retrieved. 

Thus, we can estimate that the slit flux loss for the 84 primary 
spectroscopic targets, that are perfectly centered in the slits, is 
around 15%, with a typical range between 10% and 30% de- 
pending on the seeing. For the serendipitous object with an opti- 
cal counterpart, whose off'set with respect to the center of the slit 
is known, we estimated the flux loss assuming they are compact 
under the average seeing conditions of a pointing. We show their 
distribution of the recovered flux in Figure 9. 

In order to check if the percentage of recovered flux that we 
find for the serendipitous sample is reasonable, we use a Monte 
Carlo simulation. In particular, we generate 3000 2-dimensional 
gaussian profiles with FWHM randomly extracted between 0.7 
and 1 arcsec, applying an offset between the position of the peak 
and the center of the slit randomly extracted between +1 arc- 
sec. The distribution of the recovered flux for the 3000 lines is 
shown in figure 9. We note that the choice of +1" as maximum 
offset introduces an artificial cut-off at ~10% of the recovered 




Fig. 10. Redshift distribution for the whole sample (black line), as well 
as for the different subsamples: solid lines show serendipitous lines, 
while dotted ones indicate primary spectroscopic targets; blue and red 
lines show the Ultra-Deep and deep subsamples, respectively. 



flux. This choice however is justified by the fact that an hy- 
pothetical object placed offslit by more than +1" should have 
a true Lya flux > lOx stronger than the one measured in the 
slit spectroscopy, hence its optical counterpart should be visible 
in photometry. There is an overall agreement from the result of 
this simulation and the flux loss estimated for the serendipitous 
lines with optical counterpart. For ~40% of the cases, the recov- 
ered flux is higher than 70%. The median of both histograms is 
^ 55%: we will use this value to correct the measured fluxes for 
the serendipitous lines with no optical detection in our sample. 



4. Sample properties 

4. 1 . Redshift distribution 

In Figure 10 we show the redshift distribution for our sample. 
The distribution extends from z = 2 to z = 6.7, with a median 
< z >=3. We have 5 detections at z ~ 6.5 and 9 at z ~ 5.7. The 
main VVDS targets with Lya come mainly from the Ultra-Deep 
survey, and interestingly, they show a different redshift distri- 
bution: even though the peak for the two subsample is around 
z - 2.2 2.5, the serendipitous show a broader tail at z >4, 
while 85% of the targets have z < 4. A simple Kolmogorov- 
Smirnov test rules out with more than 99% confidence that the 
two distribution are statistically equivalent. This reflects the dif- 
ferent selection for the 2 subsamples: magnitude selected for the 
former, flux selected for the latter. 

4.2. Optical counterparts and size distribution 

We searched the deep CFHTLS ugriz images, including the com- 
bined ugriz stack, for optical counterparts of the 133 serendipi- 
tous LAE. In order to account for possible positional uncertain- 
ties, we required the counterpart to be within +0.5" from the po- 
sition of the spectrum along the slit. Moreover, since even bright 
off-slit objects can show a detectable spectrum (see Sect. 3.5), 
we search for counterparts within +2" of axis. We found faint 
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Fig. 11. Full width half maximum (I'WHM) distribution as measured 
from the CFHTLS images (top panel) and from the Lyo- 2-d spectra 
(bottom panel). Targets are shown with a solid line, serendipitous ob- 
jects with and without photometric counterparts are shown with a dotted 
and dashed line, respectively. The CFHTLS image used has being built 
by stacking the 25% best-seing exposures, and has a final PSF FWHM 
of about 0.5". In the top panel, we just show serendipitous objects with 
photometric counterpart. 



counterparts for 53 (50%) and 13 (46%) LAE respectively in the 
ultradeep and in the deep sample, with magnitudes iab ranging 
from 23.5 to 27.5. These counterparts are always ±1" off axis at 
the most. For the remaining 67 objects we did not find an optical 
counterpart down to a magnitude AB ~ 28. 

In Figure 1 1 we compare the Full Width at Half Maximum 
(FWHM) for the serendipitous Lya emitters and the spectro- 
scopic targets, as measured on the deep CFHTLS images and 
directly on the 2-d spectra of the Lya line. 

In particular, we used CFHTLS-Dl (VVDS-02h field) D-25 
stack in the z-band, that has been built stacking the 25% best 
seeing images together resulting in a PSF FWHM of about 0.5 
arcsec (see Table 26 of Goranova et al. 2009). Thanks to the 
better angular resolution, these images are very useful to check 
whether or not our galaxies are resolved in the spatial direction. 
On the other side, we are also interested in measuring the size 
for our sample also on the 2-d spectra, to check the assumption 
we made in Sect. 3.5, where we estimated the slit flux loss. 

The FWHM of the serendipitous and main VVDS magnitude 
selected populations with optical counterparts measured in the 
D-25 stack are compared in the top panel of Figure 1 1 . The dis- 
tribution for both populations peaks at around 0.5 arcsec with a 
gaussian distribution width. This is expected from the measure- 



ment errors on these faint objects if they are unresolved under 
the seeing of the D-25 stack. The width of the Lya Une mea- 
sured on the 2-d spectra (bottom panel of Figure 11) indicates 
that the Lya line is emitted from a compact region, as the size 
distribution is comparable to the seeing distribution of the spec- 
troscopic observations. This is therefore indicating that most of 
the LAE are compact, both on the continuum and line emission. 

From these measurements it is clear that the faint LAE 
population is compact in size. Interestingly, there is a clear 
lack of large objects (FWHM>1-1.5 arcsec) in our sample. 
Faint LAE clearly have different sizes than the large Lya blobs 
identified in the bright LAE population (Steidel et al. 2000; 
Matsuda et al. 2004; Ouchi et al., 2008). The weak continuum 
and Lya sizes are in agreement with Venemans et al. (2005), 
Taniguchi et al. (2009), and the HST-based study of Bond et al. 
(2010). Our results support a picture where Lya emission from 
faint LAE is compact and originates from the same area as the 
UV continuum as found by Bond et al. (2010), and we do not 
confirm the trend reported by Nilsson et al. (2009) of more ex- 
tended sizes in the narrow band Lya images than in the broad- 
band images. While extended Lya emission like observed in Lya 
blobs and Lya halos is expected if Lya is produced from reso- 
nant scattering on diffuse gas surrounding the galaxy, our data 
may indicate that the gas reservoir of faint LAEs is rather com- 
pact or that, if extended, it is only scattering a small fraction of 
the Lya photons. We will explore the evolution with redshift of 
these properties in a forthcoming paper. 

We show the size distribution measured directly on the 2-d 
spectra in the bottom panel of Fig. 1 1 : due to the poorer seing of 
these observations with respect to the CFHTLS ones, this distri- 
bution is peaked at 0.9 arcsec. Hence, this justifies our choice to 
simulate lines as gaussians with an average FWHM of 0.9 arcsec 
to estimate the slit flux loss in Section 3.5. 



4.3. Spectral properties of the LAE 

We produced combined spectra for all the Lya in our sample, 
including primary VVDS targets as well as sources identified 
serendipitously. In order to preserve Une shapes and to avoid av- 
eraging out faint features, we used accurate redshift measure- 
ments from a half-gaussian fit to the red wing of the Lya lines, 
that has been used to determine the central wavelength of Lya, 
hence the redshift. Each spectrum has then been de-redshifted 
to rest-frame, normaUzed using the continuum flux in 1400- 
1800A up to z=4.6, 1400-1600A for 4.6 < z < 5.9, and 
the Lya line itself for 5.9 < z < 6.7, then averaged using a 
sigma-clipping algorithm. The results in several redshift ranges 
are shown in Figure 12 and Figure 13, for serendipitous and 
target objects, respectively. A number of spectral features can 
be readily identified at all z, including of course Lya, but also 
prominent CIV in absorption as well as Hell in emission well 
detected in the two lowest redshift bins. 

Bluer than Lya we can see that the continuum is 
more absorbed as redshift increases, with a ratio (/[1250 - 
1350A]//[1100 - 1180A] for serendipitous LAE of 1.34, 1.47, 
2.0, 9.3 for redshifts z ~ 2.3, 3, 4, 5.3 respectively, and for target 
LAE a ratio of 1 .41 and 1.51 for redshifts z ~ 2.3, 3 respectively. 
For our highest redshift objects 5.9 < z < 6.7, uncertainties on 
the continuum blueward of Lya only enables to state that this 
ratio is higher than 5. These measurements are comparable to 
the expected average extinction by the intergalactic medium, as 
Madau et al. (1995) predict ratios of 1.3, 1.5, 2.5 and 8.2 for 
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Fig. 12. Combined spectra of tiie serendipitous Lyo- emitters, in several 
redshift ranges from z=2 to z=5.9. The displayed range of flux has been 
set to show faint continuum features, while the peak of the Lya line is 
largely off-scale. Blue dotted lines indicate different spectral features 
depending on the redshift range: Ly-break at 912A , Ly/? at 1026A , 
OI at 1303A , CIl at 1333A , SilV at 1397A , CIV at 1549A , Hell at 
1640A and CIII at 1909 A . The horizontal dashed lines indicate the 
level of zero flux. 



z ~ 2.3, 3, 4, 5.3, respectively. This strengthens our identification 
of these emission line objects as LAE. 

Lyj3 - 1026 A is readily observed, and the Lyman-break at 
912A is well detected for redshifts z > 2.9 placing this feature 
in our observed wavelength range, as expected. 

To check for a difference between the targeted magnitude 
limited sample and the serendipitous LAE, we produced the 
combined spectra for all target LAE in the redshift domains 
2 < z < 2.7 and 2.7 < z < 3.5 as shown in Figure 13, 
while above z=3.5 the serendipitous LAE dominate and their 
combined spectra are quasi identical to the spectra shown in 
Figure 12. Many spectral features are clearly detected in ab- 
sorption: LyyS-1026A, OI-1303A , SUV- 1397 A . Besides Lya, 
CIV-1549A,HeII-1640A and CIII-1909A are all detected in 
emission. 

The Lyo- rest-EW increases from 38 + lA at z ^ 2.3 to 
370 + 150 at z ^ 6.3. This trend is probably the result of two 
effects: a general trend of increasing star formation at higher 
redshifts, and a difference between the observed galaxy popu- 
lations, with less (more) luminous galaxies observed at lower 
(higher) redshifts. The mean EW of Hell- 1640 for the full sam- 
ple is EW=1.4 + 0.2A , 2.0 ± O.IA and 4.0 + 2.0A at z 2.3, 
z - 3.1, and z - 4 respectively, indicating the presence of a 
young stellar population a few 10^ years old (Schaerer 2002). 
Considering only the serendipitous Lye , the EW(Hell) is ris- 
ing from 3.9 + 0.4 at z ^ 2.3 to 14.5 + 1.5 at z - 3.1. This 
is clearly indicating that a younger stellar population is present 
at the higher redshift. The difference of EW between the full 
Lya sample and the serendipitous Lya at z - 3.1 is significant, 
4.0 + 2.0 versus 14.5 + 1.5, possibly again a consequence of the 



Fig. 13. Combined spectra of the target hya emitters in the redshift 
ranges 2 < z < 2.7 and 2.7 < z < 3.5. The displayed range of flux 
has been set to show faint continuum features, while the peak of the 
Lyo- line is largely off-scale. The horizontal dashed lines indicate the 
level of zero flux. 



different populations observed, Lya in the serendipitous sample 
corresponding to fainter objects with stronger star formation or 
less extinction of the Lya emitted photons. This difference at the 
same redshift could therefore indicate that continuum-faint Lya 
contain younger stellar populations than continuum-bright Lya 
. This will be the subject of further studies. 

The Lya line in the combined spectra has an asymmetric pro- 
file with the blue wing truncated and an extended red wing, as 
shown in Figure 5. This is expected as gas outflow and inter- 
galactic absorption impose a sharp blue cutoff and a broad red 
wing. This Lya shape and the clear detection of other expected 
emission and absorption features in the combined spectra pro- 
vides further supporting evidence that the bulk of the lines in 
our sample are Lya rather than [Oil]. 

4.4. Properties of the 6 galaxies witii 5.96 < z < 6.62 

The combined spectrum of the 6 galaxies identified with 5.96 < 
z < 6.62 is shown in Figure 6. The continuum below Lya is 
barely detected in this redshift range, consistent with a strong 
absorption by the intergalactic medium, while we detect a con- 
tinuum at the 2cr level above 12 15 A. The line profile shown in 
Figure 6 is also asymmetric in this redshift range, reinforcing the 
identification of the emission lines with Lya. 

The distribution of rest EW(Lya) is in the range 70 < 
EW < 500, with a mean EW=310A and a dispersion of 150A . 
Malhotra&Rhoads 2002 suggested that EW>300Amay indicate 
a top heavy IMF. However, our measurements are just around 
this limit, and moreover very uncertain, so we cannot draw any 
conclusions on this issue. Moreover, using Eq. 4 to convert lumi- 
nosities to SFR for these 6 galaxies we obtain SFR in the range 
4 < SFR < 22MQyr^^, values that are not so extreme as to re- 
quire a top-heavy IMF. 
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Fig. 14. UGR color-color plot for galaxies in our sample between 
2.6 < z < 3.6. Filled and open circles represent respectively targets 
and serendipitous objects with optical counterpart. Black points show 
galaxies at all redshifts. The dashed box indicates the drop-out region 
in which galaxies at such redshifts are expected to lie. The cross in the 
bottom right part of the diagram shows the typical error bars. 



Fig. 15. GRI color-color plot for galaxies in our sample between 3.6 < 
z < 4.2. Filled and open circles represent respectively targets and 
serendipitous objects with optical counterpart. Black points show galax- 
ies at all redshifts. The dashed box indicates the drop-out region in 
which galaxies at such redshifts are expected to lie. The cross in the 
bottom right part of the diagram shows the typical error bars. 



None of these galaxies are photometrically detected, neither 
in the single bands reaching as faint as iab = 28 (Icr), nor in 
the combined ugirz image, reaching an equivalent magnitude of 
AB 29. 

4.5. Photometric properties 

In Figures 14 and 15 we show the ugr and gri diagrams for all 
galaxies in our sample (serendipitous and targets) at 2.6 < z < 
4.2. In fact, young star forming galaxies usually display very red 
colors around the 912A Lyman limit: the flux below 912A is vir- 
tually set to zero by the IGM absorption, while the star formation 
activity produces strong emission beyond Lya line, producing a 
very pronounced and high S/N feature. This peculiarity is ex- 
tensively used to look for high redshift (z > 2) galaxies (e.g. 
Steidel et al. 1999; Giavalisco 2002; Bouwens et al. 2007; 2009; 
2010). We determined the part of the color-color plots in which 
high-z "dropouts" are expected to be convolving synthetic spec- 
tral energy distributions for star-forming galaxies with the SDSS 
photometric system filters. 

In Fig. 14, about 80% of the galaxies with 2.6 < z < 3.6 fafl 
in the expected u-band dropout region. Similarly, in Figure 15, 
about 50% of the galaxies at 3.6 < z < 4.2 fall in the expected 
g-band dropout region. The fact that a significant fraction of ob- 
jects with secure spectroscopic redshifts falls out of the predicted 
color selection boxes has already been noted by Le Fevre et al. 
(2005b), also for LAB as Gronwall et al. (2007) find a signifi- 
cant fraction of their sample outside of a UVR LBG box, and is 
mainly the result of the photometric measurement process (Le 
Fevre et al., 2010, in prep.). 

In Figure 16 we plot, for the objects with a photomet- 
ric counterpart (primary targets and serendipitous), the equiv- 



alent width of the Lya lines as a function of the absolute 
magnitude in the UV at 1600A, for three different red- 

shift bins. Shimasaku et al. (2006), Ando et al. (2006) and 
Ouchi et al. (2008) suggested that objects with Muv < -21.5 
show typically low equivalent width (EW < 20), while objects 
with Muv > -21.5 span a wide range of EW (50 < EW < 150). 
Our data are consistent with these findings, even though some 
of our galaxies with Mf/v < -21.5 have just a lower limit to the 
EW. We note also that many of the serendipitous galaxies do not 
appear in this Figure, as they have no photometric counterparts. 

In Figure 17 we compare the Lya Luminosity, the Lya 
Equivalent Width and the star formation rate inferred from the 
Lya Luminosity for galaxies in our sample with others in lit- 
erature. The aim of this figure is to show the different selec- 
tion criteria for the various samples of Lya emitters. Narrow 
band surveys typically select objects with EWo >20 and do not 
reach Lya emitters less luminous than Liya - lO"*''^, while our 
sample has no limits on the equivalent width and goes 0.5-1 
dex deeper in luminosity. Overall, all the samples seem to fol- 
low a broad relation between Lya luminosity and Lya equiv- 
alent width. Moreover, it is clear that the typical object in our 
sample has a star formation rate SFR~ IM^yr^^ with many 
having even smaller rates, while other samples have typically 
SFR> IMovr as expected from our deeper flux limit. 

5. Luminosity functions 

5.1. Formaiism 

We aim to combine the serendipitous and target Lya emitters in 
the U-DEEP and DEEP surveys to produce luminosity functions 
at different redshifts. Since the targets and serendipitous objects 
have very different selection criteria, for the luminosity function 
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Fig. 16. Equivalent Width of the Lycf line as a function of the absolute magnitude at 1600A, for the objects with photometric counterpart,in three 
different redshift bins. Objects with detected continuum are shown as filled dots. Objects without a detected continuum, for which the equivalent 
width is only a lower limit, are shown as upperward arrows. The dotted line shows Myv = -21.5: objects brighter than this limit show a deficit of 
large equivalent widths (EWo >20). 
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Fig. 17. Luminosity versus equivalent width for the galaxies in our sam- 
ple (red points); filled and open circles indicate respectively galaxies 
with and without a detected continuum. For the latter, the measure 
of the EW is just a lower limit. Black diamonds, triangles, crosses 
and squares show respectively the samples of Ouchi et al. (2008), 
Dawson et al. (2007), Murayama et al. (2007) and Saito et al. (2008). 

calculation we treated them independently. So, each galaxy from 
the target sample has a weight calculated with respect to the 
photometric selection, and each galaxy from the serendipitous 
sample has a weight calculated with respect to the spectroscopic 
selection. 

Each of the two sample, however, is the combination of 
subsamples having different magnitude or Lya flux limits and 
covering different areas: serendipitous LAE have been identi- 
fied because their line flux exceeds the spectroscopic flux limit; 
those coming from the Ultra-Deep have a flux limit ~ 1.5 x 
lO^^^erg/ctrp-/, while those coming from the Deep survey have 
a flux limit of ~ 5 x lO^^^e rg / cm^ /s. The VVDS primary spec- 



troscopic targets are selected according to their magnitude in the 
/-band; those coming from the Ultra-Deep are selected to have 
nij < 24.75, and those coming from the Deep to have nij < 24. 

To combine together subsamples reaching different flux 
(or magnitude) limits in different areas, we used the ex- 
tended version of the l/Vmax formalism developed by 
Avni&Bahcall (1980). For each object we computed two effec- 
tive volumes Vmaxij), one for the ultradeep limit and the other 
for the deep one (if the object is a target, the limit is the mag- 
nitude limit, if the object is a serendipitous one, the limit is the 
Lyof flux limit). So, for an object with redshift z assigned to the 
redshift bin zi < z < zi- 

V„,aAj)-0(J) J —dz (5) 

where Zsup is the minimum between zi and the redshift at which 
the object could have been observed within the limits of the jth 
selection, is the solid angle covered by the jth survey, and 
dV/dz is the comoving volume element. For each object, Vmax 
is defined as: 

V,„ax = ^ V,nax(j)- (6) 
J 

This basically means that the brightest objects, that are visible 
in both the ultradeep and deep surveys, are weighted according 
to the effective volume accessible to them in both surveys, while 
the faintest one, that are visible just in the ultradeep survey are 
weighted just with respect to this latter. The target sampling rate 
and the success rate have been taken into account in comput- 
ing the volumes for the targets. They are ~ 0.1 and ~ 0.3 for 
Ultradeep and Deep surveys, respectively. 

The Lyo' luminosities have been here corrected for the flux 
slit loss, as explained in Sect. 3.5: for the targets, the expected 
slit loss is 15%, and for the serendipitous sample is ~45%. We 
verified that the choice of an average value for such coiTection 
do not affect our main results: in particular, the Schechter best-fit 
values presented in Section 5.2 do not change, within the errors, 
if the slit flux loss correction are randomly extracted for each ob- 
ject from the distribution shown in Fig. 9. Finally, we combined 
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Fig. 18. LAE Luminosity function in different redsliift bins. No IGM 
absorption correction has been applied here. Upper left, upper right 
and bottom left panels indicate 2 < z < 3.2, 3.2 < z < 4.55 and 
4.55 < z< 6.6 respectively. The error bars reflect the poissonian errors. 
Black datapoints show literature data, spanning from z=3 to z=5.7: dia- 
monds, triangles and squares represent Ouchi et al. (2008) at respec- 
tively z = 3.1, z = 3.7 and z = 5.7; asterisks, plus signs, crosses 
and circles indicate respectively Shimasaku et al. (2006) at z = 5.7, 
Murayama et al.(2007) at z = 5.7, Gronwall et al. (2007) at z = 3.1 and 
Dawson et al. (2007) at z = 4.5. We group together our measurements in 
the different redshift bins in the bottom right panel, to show that within 
the error bars, the Lya LF does not evolve from z = 2.5toz = 6. 



the two y„,ot for them using Equation 6. Moreover, a complete- 
ness correction has been applied to each object according to its 
flux and redshift, both for targets and for serendipitous emitters. 

Once we have V„uix for both the targets and serendipitous 
objects, we compute the galaxy number density for each A log(L) 
and Az bin as: 



1 



Alog(L)^ y, 



2^ 



(7) 



where n is the number of objects in that bin. 



5.2. Evolution of the luminosity functions 

We estimated the luminosity functions for 3 redshift bins: z - 
1.95 - 3., z = 3. - 4.55, z = 4.55 - 6.6 in order to keep enough 
galaxies per redshift bin for a reliable estimate of the LF. The 
results are shown in Figure 18. 

When comparing our datapoints in these different redshift in- 
tervals, we can see that the observed apparent luminosity func- 
tion (i.e. the non-IGM corrected LF) does not evolve in the red- 
shift range z - 2.5 to z ^ 6 within our error bars. Although our 
data cover a slightly wider redshift interval they are in agreement 
with the literature data, within the error bars. This is quantified 
below. 

The unprecedented depth of our survey enables to extend the 
luminosity function at z = 2 - 3. down to log(Lz,y„)=41.3, about 



one order of magnitude deeper than literature data at similar red- 
shifts. This allows us to strongly constrain the slope of the lumi- 
nosity function at z ~2.5. This parameter is extremely important, 
because a small change of this slope can produce a large change 
in the luminosity density. 

We assumed here that the luminosity function of Lya emit- 
ters is well represented by a Schechter law (Schechter 1976): 



^(L)dL = ^*(L/L*y exp(-L/L*)d(L/L*) 



(8) 



We obtained the best-fit functions for the 3 redshift bins us- 
ing this Schechter function. As a first try, we fitted the luminosity 
functions in the 3 redshift intervals allowing all the 3 parame- 
ters describing the Schechter function to vary. However, by do- 
ing this, the typical luminosity L* in the first two bins and the 
slope a in the high redshift bin are poorly constrained. Indeed, 
in the first two bins we do not have Lyo- galaxies brighter than 
lQ*^erg/cifP'/s, and in the last one we do not observe emis- 
sion lines fainter than Kf^ergjcnP' I s. So, we decided to use 
other datapoints in the literature to constrain L* in the first 
two bins. In particular, we averaged the L* values obtained by 
Ouchi et al. (2008) and Gronwall et al. (2007) at z = 3.1, obtain- 
ing Iog(L*) =42.7. So, in our fit procedure we fixed log(LLv_„) = 
42.7 in the first two bins. We verified that this approach is equiv- 
alent to including Ouchi et al. (2008) and Gronwall et al. (2007) 
datapoints with \og{Liya) > 43 to ours and performing the fit. 
Moreover, we constrained a in the last bin to the average of the 
value in the first two bins. The best-fit parameters are reported 
in Table 5.2. 

Importantly, we note that the luminosity functions reported 
above are not corrected for absorption by the intergalactic 
medium. The Lya flux from high-z sources is generally ab- 
sorbed by neutral hydrogen present in the IGM, that absorbes 
the blue wing of the Lya line, producing an asymmetric line 
profile (Hu et al. 2004; Shimasaku et al. 2006). We therefore 
measure IGM absorbed Lya fluxes, and not intrinsic Lya fluxes. 
Moreover, the amount of absorption is not the same at differ- 
ent redshifts: at high z the amount of intervening intergalactic 
medium is larger As a consequence, if the observed i.e apparent 
luminosity functions do not evolve from z ~ 2.5 to z ~ 6, the in- 
trinsic one positively evolves. The IGM optical depths have been 
estimated in various studies (i.e. Madau 1995, Fan et al. 2006 
and Meiksin 2006): all the authors agree that the amount of ab- 
sorption increases from ~15% at z ~3 to ~50% at z ~6 (at this 
redshift, basically all the blue side of the line is absorbed). Using 
the standard radiative transfer prescription by Fan et al. (2006) 
for the IGM optical depths, we can obtain the intrinsic luminosi- 
ties Lin, from the observed Lohs- At z ~2.5, 4.2 and 6 the ratio 
Lobs/Lin, is respectively 0.91, 0.73 and 0.52. 

In Figure 19 we show instead the intrinsic luminosity func- 
tions, corrected for the IGM absorption according to the pre- 
scription of Fan et al. (2006), together with our Schechter fits to 
the data. 

The first main result from this analysis is that there does not 
seem to be a strong evolution of the apparent luminosity func- 
tion between z ~ 6 and z ~ 2.5, within our error bars. In fact, 
looking at Fig. 18 we see that within our Icr errors, the lumi- 
nosity functions in the different redshift bins overlap; moreover, 
looking at the Schechter best fit parameters in Table 5.2, we can 
see again that within Icr the Schechter parameters do not evolve. 
On the other hand, by looking at the intrinsic luminosity func- 
tions in Fig. 19, we do see an evolution between z ~ 6 and z ~ 4, 
while no sizeable evolution is observed between z ~ 4 and z ~ 2. 
The observed evolution can be parametrized with the evolution 
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Table 2. Schechter Function Parameters. log(L*„^^^^) and log(L*„„) indicate respectively the value uncorrected and corrected for the IGM absorp- 
tion. Error bars give the Icr confidence levels. Values have been fixed during the fitting procedure when no error bars are shown. 



Intrinsic LF 




log(L(Lya))[erg/s] 

Fig. 19. We report here our estimates of the intrinsic luminosity func- 
tions, corrected for the IGM absorption (shaded regions). Blue, green 
and red indicate respectively 1.95 < z < 3, 3 < z < 4.55 and 
4.55 < z < 6.6 ranges. Solid lines show our best fits to the data. 

in L*, that is at z ~ 6 about 1.8 times higher than in the first bin. 
However, its significance level is only 1.5cr. 

The second interesting result comes from our ability to con- 
strain the faint end of the luminosity function. We find slopes 
a = -1.6!J}[^ at z -2.5 and -1.78+5J [° at andz ~4. Our data for- 
mally exclude a flat slope cc ~ -1 at 5 and 6.5cr, at these two red- 
shifts. Our slope values are significantly better constrainted than 
the best value of 1.36 and the marginalized value of -1.49;^q 34 
found at z ~ 3.1 by Gronwall et al. (2007), a consequence of the 
10 times fainter flux limit of our sample. Note that other authors, 
who again do not have deep enough data, do not try to fit a, 
but they rather fix it to some plausible value (-1, -1.5 and -2 for 
Ouchi et al. 2008; -1.6 for Dawson et al. 2007; -1.2 and -1.6 for 
Lemaux et al. 2009). Our analysis therefore establishes the first 
reliable estimate of the faint end slope of the luminosity function 
of Lya galaxies at z < 5. 

It is also important to state that the possible "non-Lya" emit- 
ters in our sample, described in Sect. 3.2, do not strongly affect 
these results. First of all, the luminosity distribution of the "am- 
biguous" lines (the 49 lines that cannot be unambiguously identi- 
fied as Lya based on the line ratios) is similar to that of the global 
sample, and thus it is not expected to affect the slope determina- 
tion. However, it is possible that in the 4.55 < z < 6.6 redshift 



bin the contamination is higher than in the others. In fact, for all 
the lines at A > 6920 (corresponding to z(Lya) > 4.7), we can- 
not check for other lines in the spectra, if the line is identified as 
[Oil] (see Sect.3.2). However, we showed that the contamination 
is not expected to be higher than 10%, based on the EW distri- 
bution. Thus, assuming that all the 17 lines that can be identified 
as [Oil] are at z > 4.7 would decrease the luminosity function in 
the last bin of a factor of 50% at the most. This would imply a 
weaker evolution of the luminosity function from z = 6 to z = 4. 

5.3. Evolution of the star formation rate density 

With these new constraints on the evolution of the Lya lumi- 
nosity function at these redshifts, it is interesting to estimate 
the contribution of the Lya emitters to the global star formation 
rate density of the Universe. This is not trivial, as Lya emission 
produced not only by star formation activity, but also by other 
processes like cooling radiation, AGN activity or shock winds. 
Moreover, Lya emission is attenuated by IGM and dust. 

We are computing here the SFRD using only the intrinsic 
Lya luminosity functions. In order to estimate the contribution 
of the galaxies in our sample, we integrated the luminosity func- 
tion from Liya - 0.04 xL* to log(LLy„)=44 (roughly the interval 
covered by our data). We then converted these Lya luminosity 
densities in star formation rate densities by using Equation 4. 

The derived star formation rates are shown in Figure 20, to- 
gether with the most recent estimates of the SFRD at redshift be- 
tween z=0 and z=6. Ouchi et al (2008) and Gronwall et al. (2007) 
use a sample of narrow band selected Lya emitters, while the 
other estimates are based on the UV galaxy luminosity function. 

We can see that our measurements show a slight evolution 
of the Lyof SFRD between z ~2.5 and z ~4, and a much more 
significant evolution from z ~ 4 and z ~6. Our values are com- 
patible within the errors to those obtained by Ouchi et al. (2008). 

When comparing our estimates of the Lyo' SFRD with the 
UV-derived ones, we see that the total contribution of Lya galax- 
ies to the global star formation density at z ~ 2 - 6 is important, 
increasing from ~ 20% at z ~2.5 to ~ 100% at z ~ 6. This 
implies that the Lya emission is a good tracer of the star forma- 
tion, this statistically implies that, while at redshift z ~ 2 only a 
small fraction of the galaxies contributing to the star formation 
history of the universe also show a Lya emission, all the galaxies 
at z ~ 6 do show Lya in emission. In other words, the so called 
escape fraction, that is the fraction of Lya emission produced 
by the star formation that actually escapes the star formation re- 
gions changes a lot with the cosmic epoch, from 20% at z = 2.5, 
reaching up to 100% at z ~ 6. This seems to indicate that the 
mechanism which is absorbing the Lya photons in most of the 
galaxies at z = 2 is not effective at z ~ 6, as is expected in a very 
low dust medium. However, dust estimates at z - 5 - 6 show that 
the dust content should be sufficient to produce significant Lya 
photons absorption (e.g. Bouwens et al. 2009). 
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Fig. 20. Evolution of the star formation rate density as a function of 
the redshift, inferred integrating the intrinsic Lya luminosity func- 
tions down to 0.04 x L*. Red open circles are our data, while 
blue and green open circles represent other Lyo- SFRD respectively 
by Ouchi et al. (2008) and Gronwall et al. (2007). , Open trian- 
gles, open squares, open lozenges, stars and crosses are UV esti- 
mates respectively by Schiminovich et al. (2005), Reddy et al. (2008), 
Bouwens et al. (2007), Beckwith et al. (2006), Tresse et al. (2007) and 
Bouwens et al. (2010) 



Alternatively, it is possible that the UV luminosity func- 
tion based on Lyman break galaxies searches is based on more 
and more incomplete counts at increasingly high redshifts, and 
that the current UV-derived SFRD are underestimates (e.g. Le 
Fevre et al. 2005b; Paltani et al. 2007). Current estimates of the 
luminosity density at z ^ 5 - 6 agree within a factor 2-3 or so 
(Bouwens et al., 2007, 2009), and will need new generation sur- 
veys to be improved. 

6. Summary 

In this paper we have reported the discovery of 217 faint LAE 
in the range 2 < z < 6.62 from targeted and serendipitous 
very deep observations using the VIMOS multi-slit spectrograph 
on the VLT. Adding together the areas covered by each slit- 
let combined to the wide wavelength coverage 5500-9350A of 
the Deep survey and 3600-9350A for the UltraDeep, we sur- 
veyed effective sky areas of 22.2 arcmirp- and 3.3 arcmin^ re- 
spectively. This produces a survey volume of ~ 2.5 x 10^ Mpc^, 
observed to unprecedented depth F~ 1.5 x IQ^^'^erg/s/cm^. 
This volume is about one order of magnitude bigger than all 
other spectroscopic surveys produced up to now at compara- 
ble fluxes: van Breukelen, Jarvis & Venemans (2005) sampled 
lO^Mpc^ down to lAxlQ-^'^erg/cm^/s; Martin et al. (2008) 
sampled 4.5 x 10'^ Mpc^ down to similar fluxes in a naiTow red- 
shift range. Narrow band imaging surveys sampled bigger vol- 
umes than ours, but to a shallower flux: Ouchi et al. (2008) 
covered a volume of ~ IQ^Mpc^, down to fluxes ~ 2 x 
erglcrvp- 1 s. Serendipitous surveys have been presenting 
low number of objects (Sawicki et al., 2008), even if somewhat 



deeper (Malhotra et al., 2005; Rauch et al., 2008). We are there- 
fore sampling deeper into the LAE luminosity function as we 
discussed in section 5.2. 

From an observational point of view, we demonstrate the ef- 
ficiency of blind LAE searches with efficient multi-slit spectro- 
graphs. The success of our approach is the result of combin- 
ing a broad wavelength coverage to a large effective sky area, 
with long integration times, made possible by the high multiplex 
of the VIMOS instrument. The broad wavelength coverage has 
been essential to secure the spectroscopic redshifts from one sin- 
gle observation, without the need for follow-up to confirm the 
Lye nature of the emission lines detected. This observing effi- 
ciency compares favorably with the time needed to perform nar- 
row band imaging searches followed by multi-slit spectroscopy, 
and comparing the wide range in redshift covered by the for- 
mer versus a narrow range for the latter When the density of 
faint LAE is high, of the order several LAE/arcmin^, multi-slit 
spectrographs become more efficient to secure a large number 
of confirmed sources than narrow band imaging searches, while 
at bright fluxes covering a wide field is essential to find rarer 
sources and narrow band imaging is more efficient. The two ap- 
proaches will therefore remain complementary. 

Our main findings are the following: 

1. We found a total of 217 LAE with confirmed spectroscopic 
redshifts in the range 2 < z < 6.62, 133 coming from the 
serendipitous discovery in the multi-object spectrograph slits 
of the VVDS (105 from the Ultra-Deep and 28 from the 
Deep), and 84 coming from targeted VVDS observations of 
galaxies with 17.5 < Iab < 24.75 (Le Fevre et al., 2010, 
in prep.). About 50% of the Ultra-Deep and 40% of the 
Deep serendipitous targets have a detected optical counter- 
part down to magnitude AB ~ 28 in deep CFHTLS images. 

2. The observed projected density of LAE with a Lya emis- 
sion brighter than F~ 1.5 x IQ^^^ erg I s I cnP' in the range 
2 < z < 6.6 is 33 LAE per arcmin^, with 25 LAE per 
arcmin^ with 2 < z < 4.5 and 8 LAE per arcmin^ with 
4.5 < z < 6.62. The corresponding volume density of faint 
LAE with L{Lya) > 10'*'ergs.s"' is ~ 4 x 10"^Mpc"^ a high 
density not yet observed at these redshifts. 

3. The mean rest-frame EW(LyQf ) of LAE in our sample range 
from about 40A at z ~ 2 - 3 to ~ 300 - 400 at z~ 5 - 6, and 
the star formation rate covers a wide range 0.1 - 2OM0yr"\ 
assuming no ISM or IGM extinction and a Salpeter IMF. The 
HeII-1640A emission has EW~ 4-14atz~2-4 indicat- 
ing the presence of young stars of a few Myr We therefore 
detected vigorously star forming galaxies as well as galaxies 
with star formation comparable to dwarf starburst galaxies at 
low redshifts. 

4. The Lya apparent luminosity function does not evolve be- 
tween z-2 and z-6, within the error bars of our survey. 
Taking into account the average differential evolution in the 
IGM absorption with redshift therefore translates into a posi- 
tive evolution of the intrinsic Lye LF of about 0.5 magnitude 
from z~2-3toz~5-6. 

5. We obtain a robust estimate of the faint end slope of the LAE 
luminosity function from a large sample of spectroscopically 
confirmed LAE. It is very steep: we find a - -1.6 at z ~2.5 
and 1.8 at z ~4. 

6. The SFRD contributed by Lya galaxies is increasing from 
5 - lO^^Moyr-'Mpc-i atz - 2.5 to ^ 2xlO"2M0yr-iMpc-i 
at z - 6. The contribution of the Lya galaxies to the total 
SFRD of the universe as inferred by UV luminosity func- 
tions reported in the literature increases from ~ 20% at z=2.5 
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to ~ 100% at z=6. This seems to imply that all the galaxies 
that are forming stars at z = 6 must show Lya emission, 
therefore are in a very low dust medium. At z=2.5 80% of 
the star forming galaxies must have the Lya emission pro- 
duced by star formation blocked by some mechanism so that 
only 20% of the star forming galaxies show Lya emission. 
A direct consequence would be that the Lya escape frac- 
tion varies from 0.2 at z ^ 2.5 to about 1 at z ^ 6. This 
result would remain robust only if the total SFRD estimates 
based on U V luminosity functions using Lyman break galax- 
ies identifications are complete. 

The new VVDS measurements reported here bring a new 
important constraint to the LAE LF with a steep faint end slope 
observed at 2 < z < 4.5. This implies that LAE with star forma- 
tion rates of a few 10^'Moyr"' comparable to that of low red- 
shift dwarf star forming galaxies are dominating the LAE SFRD 
at these redshifts. While it remains to be proven that the faint 
end slope of the LF stays steep going to still higher redshifts, 
assuming the slope a ^ 1 .7 as implied by our measurements at 
2 < z < 4.5 would imply that the LAE population is becoming 
the dominant source of star formation producing ionizing pho- 
tons in the early universe z ~> 5 - 6, becoming equivalent to that 
derived from Lyman Break Galaxies searches. The steep faint 
end slope further implies that during reionisation sub-L, galax- 
ies may have played an important role in keeping the universe 
ionized. 

These results further demonstrate that efforts dedicated to 
constraining the evolution of the luminosity function of high red- 
shift LAE will remain an important tool to probe into the reion- 
isation period. 
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